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SUMMARY 

A recent approach to the numerical solution of the steady Euler equations 
Is to embed the first-order Euler system In a second-order system and then to 
recapture the original solution by Imposing additional boundary conditions. • 
Initial development of this approach and computational experimentation with It 
have been based on heuristic physical reasoning. This has led to the con- 
struction of a relaxation procedure for the solution of two-dimensional steady 
flow problems. In the present report the theoretical justification for the 
embedding approach Is addressed. It 1s proven that, with the appropriate 
choice of embedding operator and additional boundary conditions, the solution 
to the embedded system Is exactly the one to the original Euler equations. 
Hence, solving the embedded version of the Euler equations will not produce 
extraneous solutions. 


INTRODUCTION 

In the development of numerical solution procedures for the steady Euler 
equations, the common approach Is to replace the steady equations by their un- 
steady counterparts and then to seek a temporally-asymptotlc steady solution, 
either In real time ([1], [2]) or In pseudo time ([3] - [5]). Due to the dif- 
ficulties associated with the numerical solution of a direct finite difference 
representation of the steady Euler equations, relatively few departures from 
this approach are to be found In the literature. Steger and Lomax [6] devel- 
oped an Iterative procedure for solving a nonconservation form of the steady 
Euler equations for subcritical flow with small shear. Desiderl and Lomax [7] 
Investigated preconditioning procedures on the matrix system arising from the 
finite differencing of the Euler equations. Bruneau, Chattot, Lamlnie and 
Gulu-Roux [8] have used a variational approach to transform the Euler equa- 
tions Into an equivalent second-order system. Preliminary numerical results 
have been presented for two-dimensional Internal flows. Recently, Jespersen 
[9] has made significant progress toward adapting multigrid techniques to the 
solution of the Euler equations and has presented results for transonic flows 
over airfoils. 
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Johnson [10] - [12] proposed a surrogate-equation technique. In which the 
first order steady Euler equations are embedded In a second order system of 
equations. The solution of the original Euler equations Is then recaptured by 
Imposing appropriate additional boundary conditions on the embedded system. 

The advantages of such an approach are that the difficulties of solving the 
direct difference representation of the steady Euler equations can be bypassed 
and the resulting second-order embedded system can be solved by a variety of 
well-proven numerical procedures. Initial development of this approach and 
computational experimentation with It have been based on heuristic physical 
reasoning. This has led to the construction of a relaxation procedure for the 
solution of two-dimensional steady flow problems. All the numerical results 
In [10] .- [12] suggest that this Is a viable and potentially more economical 
approach than the alternative of solving the unsteady equations of motion. 

In this report the theoretical justification for the embedding approach Is 
addressed. It Is proven that, with the appropriate choice of embedding opera- 
tor and additional boundary conditions, the solution to the embedded system Is 
exactly the one to the original Euler equations. Hence, solving the embedded 
version of the Euler equations will not produce extraneous solutions. The 
following section contains the proof for the two-dimensional Euler equations 
and shows that for the Cauchy-Rlemann equations a similar result follows Im- 
mediately. Generalizations to three dimensions and to other systems of par- 
tial differential equations are mentioned subsequently. 

EMBEDDING PROCEDURE 

The steady Euler equations can be written In vector form as 

fx + gy = 0 (1) 

where x and y are Cartesian coordinates. 
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Here p, p, u, v, and E denote respectively the density, static pressure, 
velocity components In the x and y directions, and the total energy per unit 
volume. Furthermore 

E = p [e + 1/2 (u2 + v2)] 

where the specific Internal energy e Is related to the pressure and density by 
the gas law 


p = (y - l)pe 

with Y denoting the ratio of specific heats. 
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By Euler's theorem on homogeneous functions, f and g can be expressed (see, 
for example, [13], [14]) as f = Aw and g = Bw, where A and B are the Jacobian 
matrices 
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Now, Eq. (1) can be written as 
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Let L denote the differential operator 
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Then the Euler equations become 


Lw = 0 (4) 

Now, let L* be the formal adjoint operator to L defined by 

L* ' - fc * I?) (5) 

where and are the transposes fo A and B respectively. We may then 
consider the Euler equations (4) as embedded In the second-order system 

L* Lw = 0 (6) 

Let D be a bounded closed region with a piecewise smooth boundary. For 
simplicity of argument, assume that Eq. (6) Is defined In a domain containing 0. 



We now show that with an additional condition on the boundary, ao, of D, solu- 
tions of Eq. (6) are also solutions of Eq. (4). 

Theorem Let L and L* be defined as In (3) and (5) respectively. If w Is a 

solution of 

L* Lw = 0 In D 

and satisfies the additional requirement 

Lw = 0 on ao, 

then It Is also a solution of 

Lw = 0 In D. 

Proof Let <• , •> denote the Euclidean Inner product In four-dimensional 
space. It can be shown that (see Appendix A for details) for any w, 

<Lw, Lw) - <w, L*Lw> = <Aw, Lw) + ^ <Bw, Lw) (7) 
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Integrating over D and using Green's theorem, we obtain 



Lw, Lw> - <w, L*Lw» dxdy 
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Here the line Integral In the last expression of Eq. (8) Is evaluated In the 
counterclockwise direction over the closed contour 80. Now, If w satisfies 
the hypotheses of the theorem, I.e. L*Lw = 0 In D and Lw = 0 on 80, then 
Eq. (8) reduces to 



(Lw, Lw)dxdy = 0 


This Implies that 


1n 0 and hence 


(Lw, Lw> = 0 


Lw = 0 


In 0. 


Q.E.O. 


Now, consider the special case of the Cauchy-Rlemann equations 


Ux + Vy = 0 

(9) 

Vx - Uy = 0 

(10) 


Let 
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and rewrite Eqs. (9) and (10) In vector form 

fx + gy = 0 


( 11 ) 


If we choose 


w = f 
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then we have 


and 
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Eq. (11) can then be expressed as 


(Aw) + ^ (Bw) 


or 


[fc <* ) * i7 <» >] “ ■ 


Hence, 1f we again use L to denote the differential operator 

the Cauchy-Rlemann equations can also be written as 

Lw = 0 
Let 


( 12 ) 
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L* = - 




ax ® ay, 


(14) 


Then Eq. (13) can be considered as embedded In 

L*Lw = 0 

Note that a few simple matrix multiplications will reduce Eq. (15) to 

2 2 

~2 ~2 w = 0 

ax'" ay*^ 

simply the well-known fact that Eqs. (9) and (10) are em- 
bedded In the second-order system (16). ' ' 

Now, let p be the same region as defined previously. Then the Introduc- 
tion of the differential operators L and L* for the Cauchy-Rlemann equations 
suggests the following Immediate consequence of the previous theorem. 

L p . rollary If w Is a solution of Eq. (16) In 0 and If, on the boundary of D, 
(lO)^ln^r^^ (9) and (10), then It Is also a solution of Eqs. (9) and 


(15) 


(16) 


6 


Thus 1f one wishes to obtain the unique solution to a boundary value prob- 
lem of the Cauchy-Riemann equations (9) and (10), one can also solve Eq. (16) 
together with the original boundary conditions and the additional requirement 
that Eqs. (9) and (10) be satisfied on the boundary.* 


GENERALIZATIONS 

Generalization of the result discussed above to the three-dimensional 
steady Euler equations is straightforward. Suppose the equations of motion 
are expressed as 


fx + 9y + hz = 0 

This can then be written as 

[fj (« ) ♦ §7 (® ) * fc >] " ’ 0 

and hence the operator L can be Introduced 

•- = i7 ^ i7 ^ il > 

The further details are analogous to the case of two dimensions. 

The embedding concept can be used on any partial differential equations 
expressible in the form 


Lq = 0 or Lq = f 

However, we shall not pursue this idea further here. 

In Oesideri and Lomax [7], preconditioning matrices are investigated. In 
our Eq. (6), L* may be considered as a preconditioning operator. Hence, the 
embedding method is a preconditioning procedure for the continuous model, 
while Oesideri and Lomax's approach is one for the corresponding discrete 
model . 

Based on the mathematical formulation presented here, two-dimensional 
steady Euler solvers are currently being developed. A detailed description of 
this work. Including numerical results, will be presented in a forthcoming 
report. 


CONCLUSIONS 

It has been proven that, for the numerical solution of the two-dimensional 
steady Euler equations, one can solve a second-order embedded system together 


*The authors understand that Dr. T. N. Phillips of ICASE - NASA Langley Re- 
search Center has recently obtained a similar result for the non-homogeneous 
Cauchy-Riemann equations. 
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with appropriate additional boundary conditions. This provides a theoretical 
justification for the recent computational experimentation with the surrogate- 
equation technique. 

The proof presented here Is extendible to three dimensions and the embed- 
ding technique Is applicable to a wider class of partial differential equa- 
tions than the Euler equations of motion considered here. 
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APPENDIX A 


Derivation of Eg. (7) For any differentiable vector-valued functions U and V, 
we have 

= |7> v> 

and 

<i7- »>=-<“> v> 

Hence, we have 

<Lw. Lw>.<|; (Aw) * ^ (Bw), Lw> 

= (fc <*"> . Lw) * (Bw) , Lw> 

' - <*“■ i7 <Lw)> t ^ <Aw. Lw> - <Bw, (LW)> * <Bw, Lw> 

= - <W. if (Lw)> t f; <Aw, Lw> - <W. B^ ^ (Lw)> . ^ <Bw. Lw> 

= -(f 'A' fc <*“• I-") * i7 I-") 

. <W, L*Lw) * <Aw. Lw) ♦ <Bw, Lw) 

This is Eq. (7). 
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